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ABSTRACT 

The well-known shock solutions of the Korteweg-de Vries-Burgers equation are revisited, to- 
gether with their limitations in the context of plasma (astro)physical applications. Although 
available in the literature for a long time, it seems to have been forgotten in recent papers that 
such shocks are monotonic and unique, for a given plasma configuration, and cannot show oscil- 
latory or bell-shaped features. This uniqueness is contrasted to solitary wave solutions of the two 
parent equations (Korteweg-de Vries and Burgers) , which form a family of curves parameterized 
by the excess velocity over the linear phase speed. 



Subject headings: plasmas - shock waves 

Among the paradigm nonlinear evolution equa- 
tions cropping up in various domains of physics, 
the Korteweg-de Vries-Burgers (KdVB) equation, 
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arises in physical media where nonlinearity, dis- 
persion and damping interact on slow timescales 
to produce solitary structures. More specifically, 
in plasma physics ([1]) typically obtains by reduc- 
tive perturbation analysis of a multi-fluid model, 
through the use of coordinate stretching 
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combined with expansions of the dependent vari- 
ables like 
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in addition to an appropriate scaling of the damp- 
ing coefficient, in many cases due to viscosity. 
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Here x and t are the original space and time co- 
ordinates, respectively, and ip refers to the elec- 
trostatic potential of the solitary waves. In the 
absence of damping (C = 0), the KdVB equa- 
tion (JTJ) reduces to the KdV equation, whereas in 
the absence of dispersion (B = 0), it recovers the 
Burgers equation, which bears kink-shaped mono- 
tonic shock profile solutions. All this is well known 
and has been in the literature for a long time, but 
we will have to come back to these points later. 

For a purely mathematical study of the proper- 
ties of the KdVB equation, ([lj is given and its 
coefficients A, B and C might be regarded as 
free parameters. However, the moment the KdVB 
equation is derived for a particular plasma (as- 
tro)physical configuration, the precise and often 
elaborate form of A, B and C has to be com- 
puted. Although the intermediate details need not 
concern us here, we still have to remind ourselves 
that A, B and C are functions of the plasma com- 
positional parameters, which also determine the 
linear phase velocity A, and thus cannot be cho- 
sen randomly. Moreover, in the process of deriv- 
ing ([TJ one has imposed/used that ip\ vanishes in 
the undisturbed medium, upstream of the shock 
or soliton solutions, translated as ip± — > for 
£ — > +oo. All this has important consequences 
for the discussion which follows. 
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Once this is properly kept in mind, there are 
several ways of deriving the stationary shock struc- 
ture of ([1]) , by changing to a co- moving frame with 
coordinate 

X = «(£ - Vt), (4) 

where n and V are related to the inverse width and 
the speed of the shock, respectively. Therefore, it 
is assumed that both k and V are positive. The 
shock solutions of the KdVB equation have been 
in the literature for a long time, and later red- 
eri ved by the so-called "tanh" m ethod, formalized 
bv lMalfliet fc Heremanl (jl996al ibh. 

However, we have to come back in explicit de- 
tail to the shock solution of the KdVB equation, 
in view of recent misunderstandings about its va- 
lidity and its applications, as shown below. One 
also has to remember that for all solitary waves, 
for which explicit analytical expressions have been 
obtained, amplitude, width (inversely related to k) 
and velocity V are inherently linked. Usually, fix- 
ing one of these parameters determines the others. 

Now, when looking at several papers in the re- 
cent literature (IShah fc Saeedl2009tlSaeed fc, Shahl 
2010; Pakzad 2OllaH0 Id uPakzad fc Javidanl201ll ; 
Shah fc Saeedl l201lt Ishah. Hague fc Mahmoodl 
20111 ). one see s that k = 1 is taken, whether ex- 



plicitlv stated (IShah fc Saeedll200fll feaeed fc Shahl 



le & y 


ahmoodl 


2011a 


bc]d|; 



Pakzad fc Javidanll201ll ). by using the shock so - 
lution in the form given bv lShah fc Saeedl (l2009h . 



No justification at all is given as to why one would 
be allowed to put k = 1, nor is there any discus- 
sion of the consequences. As we will see, taking 
k = 1 is not only needlessly stringent, but also 
erroneous, and in many cases one is not even able 
to verify that it holds, given the complexities in 
the expressions for A, B and C, except for specific 
numerical choice of all plasma parameters. Some 
other papers even leave k undetermined, as if it 
were a free parameter rtMahmood fc Ur-Rehman 
2O10HAkhtar fc Hussainboill ). 

When the transformation (j4j is applied to (TTJ, 
one finds 

d(f! dipi 
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One of the popular methods of finding the 
shock structure for (fTJ) is through the tanh 



me thod, and we will f o llow th e original paper 
by iMalfliet fc Hereman ( 1996a ). rather than a 
vast array of newcomers. We are forced to do 
so, to point out where the specific restrictions 
to plasma (astro)physics applications play a role 
and to correct some uses in the literature which 
have str ayed in this respect from t he original so- 
lutions ( Malfliet fc Heremanl 1996al ) already avail- 
able. Our treatment here i s more general than 
that of IMalfliet fc Heremanl (|1996al) , because in 
their paper A = 1 has been taken. While one can 
always rescale the absolute value of some of the 
coefficients in ([T]), one cannot easily do away with 
the sign, and we keep therefore A as determined 
by the plasma model under consideration. 

Using the transformation a = tanh x m © ancl 
noting that da/d\ = 1 — tanh 2 x, we obtain 
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Here one common factor k and one common 
bracket (1 — a 2 ) have already been divided out, to 
simplify the subsequent computations. 

The idea is then to look for solutions ipx as a 
finite power series in a, which in this case (and 
i n many others) will end wi th the quadratic term 
([Malfliet fc Heremanl H996al ). thus 



Pi = Po + Piot + foa 



(7) 



The reason that the power series breaks off comes 
from a balance between the highest nonlinearity 
and dispersive terms in . Given that the differ- 
ent powers of a are functionally independent, we 
get a system of algebraic equations, 

a : -Vfii + AfioPi - 2Bk 2 ^ - 2Ck/3 2 = 0(3) 
a 1 : -2V(3 2 + 2A(3 (3 2 + A(3f - WBk 2 /3 2 

+ 2Ck^i = 0, (9) 
a 2 : ZA^fo + QBk 2 ^ + 6Ck@ 2 = 0, (10) 
a 3 : 2Af 2 + 24Bk 2 (3 2 = 0, (11) 

determining the as yet unknown coefficients 0q, fii 
and /3 2 . Solve first (JTTJ) for j3 2 to find 
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and substitute this in ([TO]) . This allows now to 
obtain 



5A 



Solving next (|SJ) yields 



V 12Bk 2 



(14) 



Although all coefficients needed for ([7]) have 
now been determined, there is still one condi- 
tion to be satisfied before the scheme can work, 
namely ©. This was a pparently over l ooked 
or not deeme d impo rt ant dShah fe Saeedl 2009; 
Saeed fe Shahl l2010i: IMahmood fe Ur-Rehman 



20101 lAkhtar fc Hussainl 120111: iPakzadl l2011bT 



Shah fe Saeedl 120111: IShah. Hague fe Mahmoodl 
201 lUPakzadl2011bl). whi le others (|Pakzadl201 lallcl ldl: 
Pakzad fe Javidanl 1201 lh just copied the erro- 



neous solution, without going through the algebra. 
Working out ©, one arrives at 
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where for simplicity we have taken both B and 
C positive, as they usually are in most examples 
found in the literature. Adopting other sign con- 
ventions can easily be incorporated but would add 
nothing to the physics. Indeed, it is straightfor- 
ward to see that minus signs can be handled in the 
general solution by appropriate space and/or time 
reversals. Note in passing that k tanh[K(£— T^t)] = 
— Ktanh[— k(£ — W)], for any real k. 

At this stage it is clear how serious a restric- 
tion k = 1 is, for two separate reasons. First, 
all solitary wave characteristics show an inher- 
ent link between amplitude, width (inversely re- 
lated to k) and velocity V of the structure, and 
arbitrarily fixing one narrows the choices enor- 
mously. Second, assuming k = 1 means from 
(fTS")) that C — 10-B, a relation which usually can- 
not be obeyed by inserting some numbers in the 
rather complicated expressi ons B and C, as a 
glance at the papers involv e d (IShah fc Saeedj 20091: 
Mahmood fc Ur-Rehmanl l201ti ISaeed fc Shah 



2O10HAkhtar fc Hussairj201lHPakzadl2011allbUcH 
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ah 
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201 ll ) will immedi- 



Pakzad & Javidan 2011: 
Shah. Hague fc Mahmood 
ately reveal. Taken together, this implies that 
the resulting numeric s , graphs and discussions 
(|Shah fc Saeedl 120091 : IMahmood fc Ur-Rehman 



12010: Saeed & Shahl2010l Akhtar & Hussain 


2011; 


Pakzadl2011allblldldl:IPakzad fc Javidanl2011: 


Shah & Saeedl 



120111 : IShah. Hague fc Mahmoodl 1201 11 ) cannot be 
trusted. 

Using now (jT5J) in the coefficients ([12 ]) -([Ti |) 
shows that 
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At this stage the shock solution is 

3C 2 / , 2 N V 6C 2 

^ = 25AB tanh X ) + A ~ 25AB tanhx ' 

(17) 

Since B and C are assumed positive, it is the sign 
of A which will be determining the polarity of the 
kink solution. However, this only obeys the re- 
quirement that ipx — t for £ — > +oo provided one 
takes 

6C 2 
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= 24Bk' 



(18) 



Also this inherent aspect of the correct solution 
has b een overlooked in some of the recen t pa- 
pers (IShah fc Saeedl 120091 : IPakzadl l2011al lhOdL 
Pakzad fc Javidan 201 ll ). The second expression 
for V in (|18p clearly shows the link between width 
(through k) and velocity of the structure, and for 
right propagating structures V is taken positive, 
which therefore requires B to be positive. 
Finally, we arrive at the shock solution as 

3C 2 

<Pi = [1- tanh 2 X + 2(1 -tanh x )] , (19) 

where in \ we have to insert (|15[) and (|18p . giving 
C 
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The kink structure (1191) is unique, since for a given 
plasma configuration the compositional parame- 
ters fully determine A, B and C, and hence there 
is one and only one shock solution, the generic 
profile of which we illustrate in Fig. [T] once for a 
positive (upper panel), once for a negative (lower 
panel) polarity. This point has alrea dy been made 
before (|Mamiet fc Heremanl [l996al) , in a mathe- 
matical discussion, almost in passing, without re- 
ally stressing its consequences for detailed plasma 
(astro) physics problems. 
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Further remarks are in order here. Since (fT9|) 
can be rewritten as 



<Pi 




Fig. 1. — Typical KdVB shock profile, where the 
amplitude 3C 2 / (25AB) — 0.1 has been taken for 
the upper panel and —0.1 for the lower panel. 



3C 2 
25 AB 



[4-(l + tanh X ) 2 ] , (21) 



the kink is always monotonic, and no oscillatory 
part nor peak or bell-shaped curve may appear in 



ca- 
pers ( 


Shah & Saeedl200fi 


— ^ ^ — ~„ — » — i — 
Mahmood & Ur-Rehmanl 


2010; 


Saeed & Shah 2010 


Akhtar & Hussain 


201 ll: 


Pakzad 2011a b c dl;|Pakzad & Javidan 2011; 


Shah. Haaue & Mahmoo 


2011). There mav be physical situations where 



shocks including oscillatory trails or precursors 
are observed, but these cannot be described by 
the KdVB formalism. 

Note that when C = 0, the whole shock struc- 
ture disappears. This is a direct consequence of 
the very delicate balance needed between a soli- 
tary wave (KdV) and a shoc k wave (Burgers) to 
form t he combined solution ( Malfliet fc Heremanl 
1996al ). To see this more explicitly, substitute in 
(|T9")) 1 — tanh 2 x = sech 2 x, which is reminiscent of 
the typical KdV one-soliton solution. In addition, 
since reductive perturbation analysis requires that 
<pi be small enough to neglect higher-order effects, 
3C 2 /(25|AB|) should be rather smaller than 1. 

All this has to be contrasted to what hap- 
pens when C = and ([l} reduces to the stan- 
dard KdV equation, without dissipation through 
viscosity, or when B = and ([1]) becomes the 
Burgers equation, in the absence of dispersion. 
Furthermore, when C — the KdV sech 2 x soli- 
ton cannot be directly recovere d, contrary to 
what is claimed in the lite r ature (IShah fc Saeec 
20091: ISaeed fc Shahl [joifl IPakzadl 12011/ 



Pakzad fc Javidanl l201lTlShah fc Saeedl Eoill 
IShah. Hague fc Mahmoodll201l[ ). 

To see the differences, let us now first put C = 
0, return to (|8j)- (fTT|) and go again through the 
motions. It turns out that (3 2 is still given by (fT2|) . 
but fix = and (fT4"| is replaced here by 
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Hence, to arrive at the typical KdV soliton solu- 
tion in sech 2 £ = 1 — tanh 2 £, obeying ipx — > when 
£ — > ±oo, it is required that 



V = 4Bk a 



(23) 
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and now for each superacoustic soliton velocity V 
one finds a soliton of the form 



ipi = — — sech 
A 



1 V 



(24) 



Here B > is needed, which is usually the case, 
and the soliton polarity is given by the sign of A. 

Doing a similar exercise for the Burgers equa- 
tion, with B = 0, leads from JTD]) and (JTTJ> to 
/?2 = 0, in other words, (El) stop s at the linear 
term (jMalfiiet fe Heremanlll996af ) . Now © and 
© give that 



A) 



r 



/3i = - 



2Ck 



and the proper solution needs 
V = 2Ck. 



(25) 



(26) 



Taking again V as the free parameter, the shock 
solution is found as 



<Pi 



tanh 



(27) 



With the appropriate changes of notation, the so- 
lutions (|24|) and (l27l) can be found in the original 
discussion bv iMalfliet k Heremanl ( 1996a ). 

To conclude, we have discussed the intricacies 
of the proper derivation of the solitary shock struc- 
ture and its limitations in the context of plasma 
(astro)physical applications. Although these re- 
sults and restric tions have been in the liter ature 
for a long time (Malfl iet k Heremanl 1 1996al fblh it 
seems to have been fo rgotten in recent papers 
dShah k Saeedl 120091: iMahmood k Ur-Rehmanl 
2O10HSaeed k Shahl2010tlAkhtar k HussairfeoilT 
Pakzadl2011allbUdldl:IPakzad k .TayidanlSIOlllJ shah k 
" 20111 : IShah. Hague k Mahmoodl l201lh that a 



One sees that the solitary wave solutions of the 
two parent nonlinear equations, the KdV and the 
Burgers equations, are different in character, as 
they form one-parameter families of curves, de- 
pendent on the free choice of the excess velocity 
V above the linear phase speed A. 
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